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Abstract 

We investigate the physical properties of an integrable extension of the Hubbard 
model with a free parameter 7 related to the quantum deformation of the superal- 
gebra s/(2|2)( 2 ). The Bethe ansatz solution is used to determine the nature of the 
spin and charge excitations. The dispersion relation of the charge branch is given 
by a peculiar product between energy-momenta functions exhibiting massless and 
massive behaviors. The study of the finite-size corrections to the spectrum reveals 
us that the underlying conformal theory has central charge c = — 1 and critical ex- 
ponents depending on the parameter 7. We note that exact results at the isotropic 
point 7 = can be established without recourse to the Bethe ansatz solution. 



PACS numbers: 75.10.Jm, 71.10.Pm, 64.60.F-, 64.60.an 



1 Introduction 



The study of electron correlation effects in one-dimensional systems have by now attracted 
the attention of theorists for more than a half-century. The physical behavior of one- 
dimensional correlated electron models are expected to be drastically different from that 
of free electrons [1] . It turns out that the basic excitations have a collective character and 
non-perturbative techniques becomes essential. In this context, electronic lattice models 
solvable by Bethe ansatz have provided relevant insights into the physical properties 
of such systems [2j [3j 0]. Of particular interest are integrable extensions of Hubbard 
model derived exploring solutions of the Yang-Baxter equation with two fermionic and 
two bosonic degrees of freedom [5]. Representative examples are the models associated 
with the four dimensional representations of the s/(2|2) and gl(2\l) Lie superalgebras 
P, El [8]. We remark that generalizations of the Hubbard models based on the quantum 
deformations of such algebras P, Dill E] as well as on the central extension of s/(2|2) [13] 
have also been discussed in the literature. 

The purpose of this paper is to investigate the critical properties of an extended 
Hubbard model based on the quantum deformation of the twisted sl(2\2)^ algebra |10j . 
We recall here that this model appears to provide a lattice regularization of an interesting 
integrable (l+l)-dimensional quantum field of two coupled massive Dirac fermions |12j . 
Though the respective Bethe ansatz solution is known [11] it has not yet been explored 
to extract information about the physical properties of such lattice electronic model. 
Following [11] the model Hamiltonian can be re-written as, 



H = ^ Y [ C U Ci +L CT + k - C - [ l ~~ X ° n 'i~« ~ X a n i+1 - a ] + U Y n M 
i=l <r=± i=l 
L L 

+ V Y + n i,- n i+i,+) + Y Y [ c i+ c l,- c i+i,- c i+i,+ + h - c - 

1=1 1=1 

L L 

+ j Yl [ c l+4+i,- c i,- c i+i,+ + h - c - v Yl { Ui <+ + Ut '-} w 

1=1 1=1 

where c\ a and c i a are fermionic creation and annihilation operators with spin index 
a = ± acting on a chain of length L. The operator n i)(T = c\ a c i)(T represents the number 
of electrons with spin a on the z-th site. 

Apart from the standard kinetic hopping amplitude and the on-site Coulomb term U 
we see that Hamiltonian ([I]) contains additional interaction terms. They are the bond- 
charge hopping amplitudes X a and X a , the Coulomb interaction V among electrons at 
nearest-neighbor sites, the spin-spin exchange term J, the pair-hopping amplitude Y 
besides the chemical potential amplitude /i. Integrability constraints the couplings of the 
model on the following one-parameter manifold, 

X a = 1 + o-sin( 7 ), X a = 1 -CTsin( 7 ), — = V = J = Y = cos(7) (2) 
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where the anisotropy 7 is related to the g-deformation of s/(2|2)( 2 ) by q = exp[iy]. 

The potential /1 is in principle arbitrary since the model conserves the total number of 
electrons with spin a — ±. However, the invariance of Hamiltonian flTJ) by the superalgebra 
U q [sl{2\2) < ^ t> ] fixes a relation between /1 and 7, namely [TU| [TT) 

/i = 2 cos (7). (3) 

Considering the parameterization ([2j) and ([3]) one can relate the spectra of Hamiltonian 
Eq.(JT]) at the points 7 and 71 — 7. In fact, by performing a combination of particle-hole 
Cj,<j — ► c| CT and the parity q iCr — > (— 1) 4 Q )(T transformations one is able to find the following 
relation, 

if( 7 ) = 7) (4) 

Due to property (jlj) the analysis of the physical properties of Hamiltonian ([I]) subjected 
to the constraints f[2f3]) can be restricted to the anti- ferromagnetic interval < 7 < tt/2. 
In this work we shall argue that the low-energy behavior of this model in the regime 
0<7<7r/2is that of a conformally invariant theory with central charge c = — 1. The 
point 7 = is special since the model reduces to the supersymmetric isotropic s/(2|2) 
extended Hubbard model [6]. In this case it was argued that though the excitations 
are gapless the dispersion relations have a non-relativistic branch [T3J [15]. In fact, we 
found that for the electronic model ([TJ[3]) the speed of sound of the underlying low-lying 
excitations is proportional to sin(7) which vanishes in the 7^0 limit. 

We have organized this paper as follows. In next Section we shall explore the Bethe 
ansatz solution to determine the ground state and the nature of the excitations of the 
electronic model (jTJ|3]) . A particular characteristic is that the dispersion relation of charge 
excitations combines both the behavior of massless and massive degrees of freedom. In 
Section 3 we study that finite-size properties of the spectrum of the Hamiltonian (COlll) 
by both analytical and numerical approach. We argue that the critical properties are 
described by a critical line with central charge c = — 1. Our conclusions are summarized 
in Section 4. 



2 Thermodynamic limit 



Here we will determine the ground state and the nature of the elementary excitations of 
the electronic model of Section [TJ These properties can be investigated by exploring the 
diagonalization of Hamiltonian ([U13]) by the Bethe ansatz method. It was found that the 
corresponding spectrum is parameterized by the following nested Bethe equations [TT] . 



sinh(A i /2 - 17/2) 
sinh(A i /2 + z7/2) 



7V + 

n 



sinh(Aj - /ifc - ry) 
L sinh(A i - iA k + 27) 



and 



N + +N^ 

n 

fc=i 



j-j- sinh(/ij 



Atfc - 2ry) 



sinhQj - A fc - ij) _ 

sinh(/i i - A fc + ij) ^ sinh(/i i - ji k + 2i^) ' 



N:+N. 



N 4 



(5) 



(6) 
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where the integers N a denote the total number of electrons with spin a — ±. 

The eigenvalues E(L, 7) of Hamiltonian ([U13]) are given in terms of the variables Xj by, 



V(T \ 2sm (7j /-x 



To make further progress it is important to identify the distribution of roots {A.,-, /i^} on 
the complex plane which reproduce the low- lying energies of Hamiltonian (tUEl). This task 
is performed by first determining the particle number sectors of the low-lying eigenvalues. 
This is done by means of brute force diagonalization of the Hamiltonian for small chains 
L < 12 and a few values of the parameter 7. We then compare these eigenvalues with the 
results coming from the numerical analysis of the solutions of the Bethe ansatz equations 
(JS-El)- By performing this analysis we find that the ground state in the regime < 7 < n/2 
for L even sits in sectors iV + = L/2 ±1, AL = L/2 or iV + = L/2, AL = L/2 ± 1 and 
therefore it is four-fold degenerated. Due to the particle-hole symmetry it is sufficient 
to determine the respective pattern of the Bethe roots {Xj, n^} for the sector with the 
minimum possible number of roots. In Figure [1] we exhibit the ground state Bethe roots 
for L = 12 in sectors N+ = L/2, N_ = L/2 - 1 and N + = L/2 - 1, N_ = L/2. We 
clearly see that the roots Xj are real while fik have a fixed imaginary part at in/2. 
The first excited state is double degenerated and lies in sector N + = AL = L/2. In 
Figure [2] we show the corresponding Bethe roots {Xj,fik} for L = 12. By performing this 
analysis for the low-energy excitations we find that they can be described mostly in terms 
of real variables when the second Bethe roots /x^ is shifted by the complex number in/2. 
Considering this discussion we find convenient to introduce the following variables, 

A, = Af, n = \f> +i| (8) 

where G 3? for a = 1, 2. 

Now by substituting Eq.® in the Bethe ansatz equations (15116!) and afterwards by 
taking their logarithms we find that the resulting relations for are, 

I) = ^Qf ] - £ HXf - Af, I - 7 ), j = 1, ■ ■ • , N + + AL (9) 



fc=i 



and 



N+ N++N- 

,( 2 ) _ \(2) n„\ , o^n( 2 ) _ V" snP) _ \W 1 



-J2HX?-X£\2 1 ) + 2nQ?= *(Af-A«^- 7 ), i = 1, - ,JV + (10) 

fc=i fe=i 

where function $(A, 7) = 2 arctan [cot (7) tanh (A)]. 
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Figure 1: The groundstate roots Xj (crosses) and fij (circles) for 7 = ir/5 and L = 12 in 
sectors (a) N + = L/2 , iV_ = L/2 - 1 and (b) N+ = L/2 - 1 , N- = L/2. We note that 
the roots A, are the same for both sectors. 



The numbers Qj define the many possible logarithm branches and in general are 
integers or half-integers. Considering our previous numerical analysis we find that the 
low- lying spectrum is well described by the following sequence of Q ® numbers, 

Qf = _I[L-n+-n_-l]+j-l, j = 1, . . . , L - n + - n_ (11) 

+ j = l,...,^-n + (12) 

where n± are integers labeling the sector with N± = L/2 — n± particles with spin a — ±. 

For large L the number of roots tend towards a continuous distribution on the real 
axis whose density can be defined in terms of the counting function Z(A^) = Qj^/L by 



Qf 



L 

rij 

2 H 
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Figure 2: The first excited state roots Xj (crosses) and fij (circles) for 7 = 7r/5 and 
L = 12. Note that (b) has two roots Xj fixed at ±i7r/5. 



the expression, 



P 



(a)M(a) 



dZ(A}' 



(ah 



dX 



(a) 



a = 1,2. 



(13) 



In the thermodynamic L — > 00 limit the Bethe equations (MTUl) turn into coupled 
linear integral relations for the densities p^(A^) which can be solved by the Fourier 
transform method. The final result for the densities are, 



P (1) (A (1) ) 



sin 



(7) cosh(AW) 



/2) (A (2h 



7r[cosh(2A( 1 )) -cos(2 7 )] 
1 



27rcosh(A( 2 ))' ^ 
Now from the expressions for the density p ( - 1 \X^) and Eq.([7j) we can compute the 
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ground state energy per site 600(7) = limL^ooEo(L,j)/L. By writing the infinite volume 
limit of Eq. ([7j) in terms of its Fourier transform we find, 

, . f°° , cosh \uj (ir/2 — 7)1 sinhL>(7r — 7)] „ n , , 

eoo 7 =-4 sin 7 / d^ 1 1 ' ^ . , \ - for 0< 7 <^ 15 

J cosh[ci;7r/2J smh|a>7r| 2 

Let us consider the behavior of the low-lying excited states about the ground state. 
As usual these states are obtained from the Bethe equations (MTUl) by making alternative 
choices of numbers Q " 1 over the ground state configuration. This procedure is nowadays 
familiar to models solved by Bethe ansatz and for technical details see, for instance [3|[T6]. 
It turns out that the expressions for the energy e^(X^) and the momenta p^ (A^) , 
measured from the ground state, of a hole excitation on the a-th branch is given by 

POO 

£ («)( A W) =27rp (a) (A (a) ), p( a )(A (a) )=/ e {a \x)dx. (16) 

To compute the dispersion relation e^ a \p^) one has to eliminate the auxiliary variable 
\( a ) which connects energy and momentum. This is done by first computing the integrals 
in Eq. fll6p with the help of the roots densities (1141) . We then are able to eliminate the 
rapidity from e^ a \\^) and the final results for the dispersion relations are, 



£ W(pW) = 4cos( 7 )sin^Wsin 2 (^pj +tan 2 ( 7 ) 
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( 2 )(p( 2 )) = 2sin(7)sin(p( 2 )). (17) 



Note that the dispersion relation associated to particle number excitations eW(pW) has 
the interesting feature of being factorized in terms of two physically distinct types of 
dispersions. In fact, the first part has a massless behavior while the second one has 
a massive character with a mass term proportional to tan (7). By way of contrast the 
dispersion related to the spin branch (p^) is very similar to the spin- waves of the anti- 
ferromagnetic Heisenberg XXZ model. However, for low momenta the massless character 
prevails and both charge and spin excitations have a common slope at p^ = 0, namely 



e 



(«)(p(«)) ^ 2sin(7)p (a) , for0<7<- (18) 



7T 



2 



and therefore they travel with the same speed of sound v s — 2 sin 7. 

Let us turn our attention to the physical properties of the model at special point 7 = 0. 
In this case, the Hamiltonian ([IH) commutes also with the number of local electrons pairs 
[Tl] and it is proportional to the graded permutator, 

# (7 = 0) = £ (-l)^e^ 9 e^ +1) - L (19) 
j=i 0,6=1 
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where denotes 4x4 Weyl matrices acting on the j'-th site and the Grassmann parities 
are given by pi = 0, p 2 = 1, ps, = 1, and p^ = 0. 

The diagonalization of the Hamiltonian ([19]) by the Bethe ansatz was discussed in 
the literature since long ago [3j [5]. We remark that the respective Bethe equations do 
not follow immediately from Eqs. (l5ll6l) when 7^0 due to the peculiar pattern of the 
Bethe roots {/ij}. We find, however, that certain properties of the model at 7 = can 
be inferred without the need of using its Bethe ansatz solution. This is done by first 
investigating the pattern of the ground sate degeneracies of Hamiltonian ffl9|) by means 
of exact diagonalization up to L = 12. This study has reveled that the ground state sits 
in many different sectors whose total number of particles is either L or L ± 1. This tells 
us the ground state for a given L is 4L-fold degenerated and that its energy and low-lying 
excitations can be computed from the particular simple sectors iV + = L, iV_ = or 
iV + = 0, iV_ = L. Because these are typical ferromagnetic states the calculations are 
rather direct. Denoting by p the momentum of an excitation with spin a = — over the 
state N + = L, iV_ = one finds that the corresponding energy is, 

E{p) = -2L + 4sin 2 (|) (20) 

where for a finite L the momenta p = ^-K, K = 0, . . . , L — 1. 

From Eq. (l20|) we conclude that the ground state per site is 600(7 = 0) = — 2 and that 
for low momenta p the excitation energy are proportional to p 2 . Therefore, the system 
has a nonrelativistic behavior in accordance with previous works in the literature [T4"t [15]. 
Interesting enough, we observe that such results can also be derived from Eqs. (115111 71) by 
taking the limit 7 — > 0. To obtain the ground state energy from Eq. fll5|) we first perform 
the change of variable u — > cj/7 and afterwards take the 7-^0 limit. On the other hand, 
the dispersion relation e(p) = 4sin 2 (p/2) follows directly from Eq. ([T7|) by substituting 
7 = 0. 

We have now the basic ingredients to investigate in next section the finite-size effects 
in the spectrum of the electronic model (IIH) for < 7 < tt/2. 



3 Critical properties 

The results of previous Section suggests us that the generalized Hubbard model ([U12]) in 
the regime 0<7<7r/2is conformally invariant. This means that the corresponding 
critical properties can be evaluated investigating the eigenspectrum finite-size corrections 
[T7] . For periodic boundary conditions, the ground state E Q (L, 7) are expected to scale 
as, 

L ~ Co ° 6L 2 

where c is the central charge. 



+ 0(L- 2 ), (21) 
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From the excited states E a (L,j) we are able to determine the dimensions A a (7) of 
the respective primary operators, namely 



E a (L n ) E (L^) _ 2nv B { 1 )X a { 1 ) 
L L L 2 



+ O (L" 2 ) . 



(22) 



A first insight on the structure of the finite-size corrections can be obtained by applying 
the so-called density root method [TBI UHl [20]. This approach explores the Bethe ansatz 
solution and it makes possible to compute the O {L~ 2 ) corrections to the densities of roots 
pW(AW). This method is however only suitable for systems whose ground state and low- 
lying excitations are described by real roots. Fortunately, this is exactly the situation 
we have found in Section [2] once the second root is shifted by in/2. Considering this 
subtlety on the root density approach we find that the leading finite-size behavior of the 
eigenenergies is, 



E{Ln) , . 2vr 

— ^ — = 600(7) + jpv.o) 



1 -ym,m. 

6 n +> n - 



(7) 



+ O (L- 2 ) 



(23) 



where the dependence of the scaling dimensions A^™~ (7) on the anisotropy 7 is 



-ym,m. 
n^- ,n. 



n 2 + n 2 + 2(1 - —)n 

71 



+ 



7i 



47(7T — 7) 



m 



m 



2(1 

V 71 ' 



\mm- 



As before the integers n± parameterizes the numbers of electrons N± 

m 



(24) 

L/2 — n± with 

spin a — ±. The indices m = m + + m_ and m + characterize the presence of holes in the 
and Qj distributions and in principle can be integers or half-integers. This approach 
is however not able to predict either the possible values for the vortex numbers m and m + 
as well possible constraints with the corresponding spin- wave integers n and n + . To shed 
some light on this problem we shall first study the finite-size effects at the particular point 
7 = 7r/2. For 7 = 7r/2 we see that all the interactions in the Hamiltonian ([T][3]) cancel 
out and we remain with two coupled free fermion models. In this case standard Fourier 
technique is able to provide us the exact expressions for the low-lying energies in the case 
of arbitrary L. The respective calculations depend on the total number of electrons on 
the lattice L. We find that when n = n + + n_ is odd that the expression for the lowest 
energy in this sector is given by, 



E dd(L, 



71 . 



[cos(^t) +cos(^)] 



sin 



(i) 



(25) 



Considering the asymptotic expansion of Eq. ff25l) for large L one finds 



L 2 



vJ-, 



1 
6 



n 2 , + n 2 1 



O (L- 2 ) . 



(26) 
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By way of contrast when n 



n. 



is an even number the lowest energy is, 



whose expansion for large L is, 



-7i\ 





7r(n CT +l) 


cos 


L 





7r(n CT — 1) 


COS 


L 



sin 



2tt 



(f) 



sin 



(!) 



(27) 



2 



1 n 2 + n 2 _ 1 



+0 (L~ 2 ) 



(28) 



Taking into account Eqs. (l26ll28l) we see that the expected finite size corrections depend 
whether the index n is an odd or even integer. In addition, by comparing Eqs. (l26|l28l) 
with the general results ( I23T241 at 7 = tt/2 we clearly see that for n odd the numbers m 
and m_ appear to start from zero while for n even the lowest allowed value for m and 
m_ is in fact one-half. This analysis strongly suggests that possible values for the vortex 
numbers m and m + should satisfy the following rule 



for n odd 
for n even 



0,±1,±2, 



1 3 5 
m,m_ = ±~,±-,±~, 

2 2 2 



(29) 



Let us now check if the above proposal remains valid for other values of the parameter 
7. This is done mostly by solving numerically the original Bethe equations (JU E]) up to 
L = 32. For the excited states whose respective Bethe roots are unstable already for 
moderate values of L we have used the data obtained from the numerical diagonalization 
through the Lanczos method. This numerical work enables us to compute for each L the 
following sequence 

By extrapolating X(L) for several values of L we are able to verify the expression (1241) 
for X^^il) an d the constraints (1291) . In Tables [Tj, [H and [3] we exhibit the finite-size 
sequence ( 1301) for six lowest dimensions on the even sector to make an extensive check 
of the less unusual part of the rule fT29l) . For sake of completeness we also present three 
conformal dimensions corresponding to the n odd sector. All those numerical results 
confirm the conjecture (|2"4"1 121?j) for the finite size properties of the generalized Hubbard 
model flUED- 

We shall now proceed with a discussion of the results obtained so far. From Section 
[2] we know that the ground state sits in the sectors n+ = ±1 and n_ = or n + = and 
n_ = ±1. Considering the rule (1291) the corresponding vortex numbers have the lowest 
possible values m = m + = and from Eqs. f)23ll24p we derive the following finite size 
behavior, 

E (L,j) ^(7) 



+ 



6L 2 



+ O (L- 2 ) 



(31) 
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T 

Li 


i i 

Y 2 ' 2 / 7T \ 

A o,o UJ 
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A o,o I5J 


i i 
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Q 
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U.oloooU 


l.zz f 954 


n /i 00^70 
U.4ootuz 


U.ooUzol 


n tix on x © 
U. /OZyoo 


U.bolU/o 


1 9 




1 9^Q1 90 




^77^1 n 


D 7M 9^0 




16 


0.312787 


1.243526 


0.502997 


0.376297 


0.750687 


0.699292 


20 


0.312689 


1.245526 


0.505395 


0.375829 


0.750434 


0.701996 


24 


0.312633 


1.247012 


0.506892 


0.375575 


0.750299 


0.703591 


28 


0.312595 


1.247954 


0.507906 


0.375423 


0.750218 


0.704619 


32 


0.312576 


1.248063 


0.508633 


0.375323 


0.750166 


0.705327 


Extrap. 


0.31249±1 


1.2504±2 


0.51219±1 


0.37498±1 


0.74999±1 


0.708336±1 


Exact 


0.3125 


1.25 


0.5125 


0.375 


0.75 


0.70833. . . 



Table 1: Finite size sequences ( 1301) of the anomalous dimensions for 7 = 7r/5, it/ 3 from 

i i 

the Bethe ansatz. The expected exact conformal dimensions are X 2 q 2 (7) = 4( - 1 _ 1 7 ^ , 
i _i ^ i i ^ 

^,0 2 (7) = 4(^0' ^1,-1(7) = i + 4(1- 7 /tt) ■ 



L 




^S-i(f) 




*So°(S) 


*2°:-i(i) 


1 1 

Y 2 ' 2 I 7T \ 

A 2,-2l3' ) 


8 


0.251098 


0.642630 


1.000395 


0.252587 


0.902529 


1.541194 


12 


0.250523 


0.646574 


1.047807 


0.251149 


0.910140 


1.622978 


16 


0.250301 


0.648003 


1.068757 


0.250646 


0.912924 


1.655748 


20 


0.250195 


0.648689 


1.080196 


0.250413 


0.914244 


1.672266 


24 


0.250136 


0.649072 


1.087265 


0.250287 


0.914971 


1.681827 


28 


0.250101 


0.649307 


1.092001 


0.250211 


0.915414 


1.687896 


32 


0.250077 


0.649463 


1.095375 


0.250161 


0.915704 


1.692009 


Extrap. 


0.250003±1 


0.65003±1 


1.1124±1 


0.250004±2 


0.9167±2 


1.70825±1 


Exact 


0.25 


0.65 


1.1125 


0.25 


0.91666... 


1.70833. . . 



Table 2: Finite size sequences (1301) of the anomalous dimensions for 7 = 7r/5, it/ 3 from 
the Bethe ansatz. The exact conformal dimensions are ^i'o(t) = \, ^2-1(7) = \ + 
y\\ <„\ _ ±r 1 1 

^2,-217 J - 7T T" 4(1-^/77)- 
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i i 
*iT(l) 




4 


1.497964 


1.395242 


1.426027 


1.355168 


1.225412 


1.261369 


6 


1.661634 


1.288542 


1.741570 


1.480050 


1.172763 


1.563850 


8 


1.724314 


1.224764 


1.888707 


1.523739 


1.140033 


1.639713 


10 


1.754379 


1.188829 


1.947347 


1.544502 


1.121710 


1.679436 


12 


1.771238 


1.167328 


1.977593 


1.556018 


1.110818 


1.699900 


14 


1.781660 


1.153629 


1.995956 


1.563070 


1.103911 


1.712660 


16 


1.788561 


1.144426 


2.008012 


1.567703 


1.098013 


1.721046 


Extrap. 


1.812±1 


1.12±1 


2.06±1 


1.583±1 


1.079±1 


1.73±1 


Exact 


1.8125 


1.1125 


2.05 


1.5833... 


1.0833... 


1.75 



Table 3: Finite size sequences ( I30p of the anomalous dimensions for 7 = ir/5, 7r/4 
from Lanczos. The expected exact conformal dimensions are xI'q^) = \ + 4( 7 / 7r )(i_ 7 / 7r ) > 

XI? (l) = (1 -7/vr) + i ( T ^y, X^Hl) = (1 - -y/n) + 



Direct comparison between Eq.( l2Tj) and Eg. (1311) leads us to conclude that the central 
charge of the underlying conformal theory is, 

c = -1 for < 7 < - (32) 

The conformal dimensions of the primary operators X™'™+ (7) depend on the anisotropy 
7 and they should be measured from the ground state Eo(L,j). Considering Eqs. (123f24p 
together with Eq. (l3Tl) we find that they are given by, 

X^ + + {l)=X^{ 1 )-\ for 0<7<| (33) 

To our knowlodge, models exhibiting this kind of universality class have so far been 
found in a not self-adjoint theory based on the deformed osp(2\2) symmetry [21]. There- 
fore, the correlated electron system ([U13]) appears to be the first example of a Hermitian 
Hamiltonian whose continuum limit is described by a field theory with c = — 1 with con- 
tinuously varying anomalous dimensions. The fact that a line of critical exponents with 
c < can be realized in terms of Hermitian models could be of importance for practical 
applications in condensed matter such as in the physics of disordered systems. 



4 Conclusions 

We have studied the physical properties of an exactly solvable generalization of the Hub- 
bard model with free parameter 7 related to the quantum U q [SU{2\2) superalgebra where 
q = exp(i'~f). We have determined the nature of the ground state and the behavior of 
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the elementary excitations. The peculiar feature of the model is that the dispersion re- 
lation for the charge sector is given in terms of the product of massless and massive 
energy-momenta relations. In the regime < 7 < 71/ 2 the low- lying excitations have a 
relativistic behavior and the underlying conformal theory has central charge c = — 1 and 
a line of continuously varying exponents. For the particular point q = 1 we have argued 
that basic properties can be obtained without recourse to the Bethe ansatz solution. We 
expect that this observation remains valid for all integrable models based on the Lie su- 
peralgebra sl(p\q) for arbitrary finite number of bosonic (p) and fermionic (q) degrees of 
freedom. This suggests that the models based on the deformed sl{p\q) symmetry may 
also have excitation modes with dispersion relation exhibiting both massless and massive 
behaviors. We hope to investigate this interesting possibility as well as its consequences 
in a future publication. 
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